Entrainment and Unit Velocity: Surprises in an Accelerated Exclusion Process 
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We introduce a class of distance-dependent interactions in an accelerated exclusion process (AEP) 
inspired by the observation of transcribing RNA polymerase speeding up when "pushed" by a trailing 
one. On a ring, the AEP steady state displays a discontinuous transition, from being homogeneous 
(with augmented currents) to phase-segregated. In the latter state, the holes appear loosely bound 
and move together, much like a train. Surprisingly, the current-density relation is simply J = 1 — p, 
signifying that the "hole-train" travels with unit velocity. 

PACS numbers: 05.40.-a, 64.60.an, 05.70.Fh 



Over four decades ago, the totally asymmetric simple 
exclusion process (TASEP) was introduced by two dis- 
tinct communities: biochemistry and pure mathematics 
P~H5]. This venerable model has since enjoyed much at- 
tention, especially from statistical physicists |H [5] . Vari- 
ations to the original TASEP emerged, as different fea- 
tures are recognized to be crucial for capturing essential 
aspects in biological and physical systems. For exam- 
ple, for applications to various transport phenomena in 
molecular biology [HJ [7] , these additions include extended 
objects and inhomogeneous hopping [TJ [2J [5HTT], non- 
conserving particle numbers and multiple species |121ll3j . 
"recycling" and competition [TIHTB] . Inspired by the co- 
operative increase in speed in transcribing RNA poly- 
merases (RNAP) that "push" each other [Mil]], we set 
out to study an extension, in which the particle arriving 
at the rear of a cluster of particles "triggers" the parti- 
cle at the front to move (FigjlJ. A mechanical "push" 
example is Newton's balls. In general, the mechanism of 
"pushing" need not to be mechanical or involve actual 
collisions. 

While interacting driven lattice gases have been stud- 
ied for nearly 30 years [2TJ [22] , such "facilitated" action 
was considered only recently [123-26 . Our model differs 
substantially from these previous studies: In |23j . neigh- 
boring particles may attract /repel each other or a parti- 
cle may move as far as two sites, but do not trigger oth- 
ers at a distance to move. Implied by the original name 
- restricted asymmetric exclusion process [21], particles 
(with non-zero headway) in the "facilitated asymmetric 
exclusion" process [2S] do not move unless there is one or 
more particles behind it. In the more general "coopera- 
tive exclusion process" [2H] , these rules are used a fraction 
of the time while ordinary TASEP rules apply otherwise. 
Thus, it is natural that the average current J (as a func- 
tion of p, the particle density on a ring) is always lower 
than the standard Jtasep = p0-~ p)- By contrast, in 
our model J is always higher. Thus, we name our model 
the "accelerated exclusion process" (AEP). Beyond this 
expected increase in J, we discover a qualitatively new 
phenomenon. Typically, there are two branches of J (p): 



An augmented current (AC) branch at lower densities, 
and a branch of lower current at high densities, with a 
discontinuous jump as p crosses some critical value. Even 
more remarkable is that J = 1 — p in the latter branch, 
i.e., the system has unit velocity (UV) on the average 
- with particles and holes in opposite directions. Such 
a simple property emerges from a system with interact- 
ing particles is astonishing; it behooves us to name this 
state the "UV phase." Preliminary studies indicate that 
its presence depends on two competing factors, to be de- 
tailed below. In this letter, we first provide the biological 
motivation and define our model. Simulation results, as 
well as some analytic understanding, will follow. In a 
final paragraph, we outline the issues which should be 
pursued more systematically. 
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FIG. 1. AEP with ^ m ax = 3. (a) Backward hops and overlaps 
are forbidden, (b) A particle hopping to the back of a 3-cluster 
triggers an additional hop. No other hops are triggered, e.g., 
(c): particle hits a cluster with £ > ^ ma x- 

Motivations and Model: Our model is motivated by 
a cooperative effect in transcription where the forward 
motion of an RNAP is accelerated by the presence of a 
trailing one that prevents the first RNAP from entering 
alternative kinetic pathways such as pausing and back- 
tracking 19, 20 . Although the pausing and backtracking 
aspects were incorporated in an exclusion process earlier 
|27j . facilitated motion was not studied [25]. Another in- 
stance of pushing in transport is tailgating in ve- hicular 
and pedestrian traffic, where the acceleration is mediated 
by the transfer of information (rather than mo- mentum 
as in Newtons balls). Here we model the generic coop- 
erative aspect of pushing with the following rules: Con- 
sider a 1-D lattice with periodic boundary conditions - 
for simplicity and as a baseline study. Let the occupa- 
tion variable of site i (= 1,...,L) be ni — 1 or 0, cor- 
responding to a particle or a hole. In standard TASEP, 
the system evolves by discrete attempt steps: Picking 
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a random occupied site and moving the particle to the 
next site provided the target is empty. The total number 
of particles (N = E,rij) is a constant, and the particle 
density (p = N/L) acts as a control parameter. When 
the system settles down into a steady state, one of the 
non-trivial quantities of interest is the average current, 
J (average number of particle hops in an attempt) and 
its dependence on p. Clearly, Jtasep = (n% (1 — nj+i)), 
where the average is over all possible configurations {rii} 
with the appropriate weight. Since Spitzer 3J showed 
that all {rii} (with a fixed N) in the steady state are 
equally probable, finding an exact expression for J (p; L) 
is simple [5] and Jtasep = P (1 — p) emerges in the ther- 
modynamic limit. 

For AEP, we include an additional move to account 
for long-range interactions: If a particle moves to a va- 
cant site and becomes the rear of a particle cluster, 
it "triggers" the particle at the front of the cluster to 
move in the same attempt. Note that the second par- 
ticle does not trigger another move, even if it arrives 
at the back of another cluster. Since particles are in- 
distinguishable, we can also regard this action as the 
chosen particle being "accelerated" through the entire 
cluster it hits! An exact expression for the current is 
Jaep = (th (1 - n i+ i)) + (rii (1 - n i+2 ). Being typ- 
ically higher than Jtasep > the term "augmented current" 
(AC) seems apt. It is reasonable to consider finite range 
for such interaction, an aspect we implement by allow- 
ing only clusters of length up to £ max to facilitate the 
accelerated move summarized in FigJT] This length scale 
is dependent of the specific system at hand and more 
realistic rules can be introduced, e.g., the rate of acceler- 
ation being a smooth function of £ rather than the step 
function. Although Jaep (ft) > Jtasep (p) is still ex- 
pected, simulations reveal the existence of AC and UV 
branches, with a discontinuous transition when a critical 
density, p c (^max), is crossed (Fig{2]). Before presenting 
these data, we first comment on two other perspectives 
of the model. 

Unlike TASEP, AEP is not particle-hole symmetric. 
The holes' perspective proves rather useful. When a 
hole at site i is chosen and can move to site i — 1, it 
"pulls along" the next hole, provided the latter lies within 
[i + 2, i + 1 + £ max \. A picturesque way to regard this 
action is that of a train, a language to which we will re- 
turn extensively. A third alternative description is the 
zero range process (ZRP) [2§], to which TASEP can be 
mapped. On a 1-D ring of sites labeled by a = 1,...,H, 
each can be occupied by an unlimited number of parti- 
cles, £ a , piling up in a column. A random site is chosen 
and the particle at the top is moved to the next site. It is 
clear that the sites in ZRP represent the holes in TASEP 
(with H = L — N) , while £ a represents the cluster of par- 
ticles behind hole a. The modification for AEP is simple: 
A particle at site a moves by an extra step (i.e., to a + 2) 
provided £ a+1 € [Mraax]- 

Simulations and Analytic Understanding: Deferring a 
more systematic investigation of AEP to another publica- 



tion [SU], we report the highlights of our findings, mainly 
for a L = 1000 ring with a range of p and £ max . The 
average current J is measured by the total moved parti- 
cles per Monte Carlo step(MCS) over 10 6 MCS. Figga) 
shows that Jaep (p) is indeed > p (1 — p) everywhere 
and displays more complex behavior. For £ max ^$ 10 
or £ max = L j it is a relatively smooth function, much 
like Jtasep- However, for 10 < £ max < 900, there is a 
sharp transition to a lower current branch. More remark- 
ably, this regime is well described by J = 1 — p i.e., the 
holes hop as if they are non-interacting at unit velocity! 
In summary, Jaep follows one of two £ max -independent 
functions: Jac (p) or Juv = 1 — P- The jump from one 
to the other is quite sharp. Located at some p c (£max), 
it appears to approach a discontinuity singularity in the 
thermodynamic limit. As the fuzzy lines in Fig[2] im- 
ply, large fluctuations are associated with the transition 
region, the details of which remain to be systematically 
studied. 




FIG. 2. (a) Jaep(p) for various Umax's. Detailed views of 
^max = 20, 100, 500: symbols in (b), with dashed lines for Jac 
and Juv and in (c), J(p) + p is plotted to accentuate the 
transition into UV. 

AC Branch: Given that two particles can move in the 
same attempt, the roughest estimate for Jaep would be 
2 Jtasep- Indeed, Jac is qualitatively so, yet subtly dif- 
ferent, as shown in Fig[2](a). Similar to the A < 1/2 sys- 
tems in [26] , Jac (p) also has a region of truly accelerated 
or "facilitated" motion, i.e., d 2 J/dp 2 > 0. Unlike those 
systems, Jac is much larger, while the facilitated region 
is limited to p < 0.2. Other notable features of Jac in- 
clude: Its maximum occurs at a density beyond 1/2; it 
exceeds 2 Jtasep for most of the p > 1/2 region; and the 
inflection point appears to coincide with the maximum 
of 2 Jtasep — Jac- Understanding these features remains 
a challenge [3D] . 

"Hole-train" in UV: To gain some insight on the strik- 
ing UV behavior, we begin with a simpler system: An 
infinite lattice completely filled except for H holes and 
4nax ~* 00 ■ After presenting the exact solutions for 
H = 1,2,3, we provide a general picture of "entrain- 
ment." In this scenario, on average, the holes are loosely 
bound to some finite length and move together with UV. 
This picture naturally evokes the term "a hole-train." 

Clearly, a single hole moves with UV. With two holes, 
let £ denote the gap between them. If the first hole (on 
the left in our model) is chosen for update, only the first 
moves if £ = 0, while both holes move if £ > 0. If the 
second is chosen, the first stays and £ decreases by unity 
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(if £ > 0). Since a positive £ can never increase, the 
system ends up in an "absorbing state" with £ — or 
1. In other words, the two holes form a tightly bound 
state, with equal probability to be in either £, while the 
average velocity is easily found to be unity [30] ■ We will 
refer to such a pair as an "engine" (of a train) . The first 
non-trivial case is H — 3. Since the third hole has no 
effect on the first pair, an engine will eventually form 
and proceed with UV. The third hole, naturally named 
the "caboose," can trail the engine by any distance, to. 
Thus, a configuration is uniquely specified by £ (0 or 1, as 
above) and to, and a complete solution is given in terms 
of the probabilities p{£,m). Writing a master equation 
for p and following standard generating function tech- 
niques, we find that, apart from the first few terms, the 
stationary p's decays as C™\ with C = 2 — \/2, so that 
(to) ~ £/ (1 — () = v2 [30]. In other words, the caboose 
is also bound, though not as tightly as the first two holes. 
It trails with an exponential tail of characteristic spacing 
/i = — 1/lnC — 1-85. We can show rigorously that this 
3-hole-train moves with UV on average. Proceeding to 
H > 3, we can provide the following convincing argument 
(as an exact solution is being sought) for entrainment. 
The total train length - from the first hole to the last - 
can change only when three particular holes are chosen 
for update: The first, the last, and the next-to-the last. 
Choosing the first always increases the length by unity. 
If one of the other two is chosen, the length may either 
remain unchanged or decrease by unity. Naively then, for 
every chance to lengthen the train, there are two which 
may shorten it! This length cannot fall below H and so, 
the expected exponential tail also emerges! Moreover, as 
the engine necessarily moves with UV, the whole train 
must also move at UV. In the ZRP framework, this pic- 
ture is even easier to grasp. The system corresponds to an 
open lattice of H sites, with unit entry rate - regardless 
of the contents of the first site. At the opposite end, the 
exit rate may be as high as two - when the last two sites 
are both occupied and chosen for update. If the filling 
fraction happens to be high, the imbalance of entry-exit 
rates will bring it down. Thus, the overall density is ex- 
pected to remain finite. Meanwhile, being limited by the 
entry rate, the average current is precisely unity. 

Next, we turn to applying these results to our ring 
where £ max ,L < oo. While UV for H = 1 is trivial, the 
non-trivial role of £ max already emerges when H = 2. If 
f max > L, every attempt results in both holes moving 
(except for I = 0, which quickly becomes i — 1). Indeed, 
I never changes and we have 4 moves in one MCS, re- 
gardless of which particle is chosen. Thus, we have many 
"absorbing states," each with current (4/L) that even ex- 
ceeds 2J T asep = 2 (2/L) (1 — 1/ (L— 1))! On the other 
hand, for small £ max (say, 10), two holes far apart will 
perform independent (totally biased) random walks so 
that there are 2 moves per MCS. However, fluctuations 
will cause the smaller of the two gaps between them to 
fall below ^ max . From this point on, i cannot increase, 
as the lead hole will always pull the trailing one. This 



scenario continues until the pair forms an engine moving 
with UV. Here, the finiteness of L is irrelevant. 

For H = 3, the role and value of £ max becomes more 
significant, since the caboose is loosely bound and /i must 
enter somewhere. In particular, if £ max is too small, the 
caboose can easily come unbound, wanders around the 
ring and then "unbinds" the engine by pulling the first 
hole away from the second. On the other hand, if £ max is 
too large (e.g., O(L)), the caboose directly affects the 
engine's integrity. A clear picture now emerges: For 
/i <C £max <C L, the 3 holes tend to be entrained and move 
with UV. Indeed, simulations with ^ max = 20 shows that, 
in 4000 measurements, the train length never exceeds 15, 
a fact entirely consistent with £ 15 ~ 10~ 4 . This picture 
extends easily to 3 < H <C L as a hole-train moving at 
UV for a range of £ max . To provide a more quantitative 
view of the role played by £ max , we show various clus- 
ter size distributions in the H = 5 case in Figj3j With 
£max — (ordinary TASEP), we find a broad distribution 
of sizes, implying the absence of entrainment. The aver- 
age cluster size, as expected, is ~ 200 = L/H. At the 
other extreme (£ maK = 1000), we find a statistically indis- 
tinguishable distribution! This result is also understand- 
able, especially from a ZRP perspective, where particles 
just move around the five sites, at typically twice the 
speed. The other three distributions are drastically dif- 
ferent, showing dominant peaks at both ends, a marked 
signal of strong clustering of the 5 holes. For £ max = 10 
and 20, the system is clearly attempting a transition to 
the entrained state. There is a small but broad distribu- 
tion of sizes in between, implying that one of the holes 
becomes unbound, creating intermediate size gaps as it 
wanders around the rest of the ring. By £ max = 100, such 
events do not occur in our runs. Indeed, the frequency of 
the large cluster shows that just one such cluster appears 
in each measurement. Of course, its size is precisely the 
complement of the length of the hole-train. In this case, 
the latter is ~ 8.92, approximately 2H — 1. In short, the 
hole-train has, on average, one spacing between its cars. 
We will show in a more systematic study that this result 
prevails in the UV phase for a large range of H [501. 




10" 
cluster size 

FIG. 3. Cluster-size distribution for H = 5 in L = 1000. 
Broad distributions for the ordinary TASEP and AC are es- 
sentially identical. In phase-segregated UV, dominant max- 
ima prevail at both extremes. For quantitative comparisons, 
the TASEP distribution (black line) is shown in all panels. 

Given the existence of a UV phase, let us turn briefly 
to the transition to the AC phase. From the viewpoint of 
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the hole-train, it is an "unbinding" transition. As more 
holes are added to a finite ring, the size of the macro- 
scopic cluster separating the engine from the caboose, A, 
decreases. When it drops to ~ £ max , the caboose will 
"pull the engine apart." Meanwhile, if we assume the 
train length to be ~ 2H, we have A « L — 2H, leading to 
a rough estimate p c ~ (1 + Cax/^) /2. Remarkably, this 
estimate is within 10% of the data points gathered so far 
[50] . In the transition region, the train dissolves and re- 
forms, resulting in the large fluctuations we observe (in 
Figj2|. To find a good estimate for p c is non-trivial since 
these fluctuations will undoubtedly play important roles. 

Concluding Remarks: Inspired by long-range interac- 
tions among particles such as the speed-up in transcrip- 
tion through cooperative RNAPs, we investigate AEP in 
which a particle hopping onto a cluster of length up to 
f max simultaneously triggers the first particle in that clus- 
ter to hop. This extension from the paradigmatic TASEP 
gives rise to various novel properties such as the tran- 
sition from homogeneous to phase-segregated and the 
intriguing unit-velocity phase. Simulating such a sys- 
tem on a L = 1000 ring with various filling fractions 
p, we find the augmentation of J (p) over the ordinary 
>/tasep = p (1 — p)- Surprisingly, we discover that when 
p exceeds a critical p c (£ max ), a "condensation" transi- 
tion takes place and the system becomes phase segre- 
gated. Here the holes gather into a loosely bound clus- 
ter, moving as a whole around the ring, motivating us to 
name it a "hole-train." Even more remarkably, this train 
moves with unit velocity! Focusing on the UV phase, 
we measured cluster-size distributions and sought theo- 
retical understanding. Our studies with small number 
of holes provided adequate insight for us to understand 
why entrainment exists, why the train should move at 
UV, and how the transition from the homogeneous AC 
phase arises. The role of £ max is critical: If it is too 
small, binding cannot be sustained; too large, the two 
ends of the train interact and the engine/train disinte- 
grates. Perhaps these insights will help us arrive at a full 
analytic theory. 

Many other intriguing issues remain, on both theoreti- 
cal and modeling fronts. We should first emphasize that, 
at this stage, AEP should be viewed more as a significant 
extension of TASEP than an explicit model for transcrip- 
tion. To address the former, the most immediate need 



is a systematic study with a range of L's so that finite 
size effects can be quantified. Locating the transition 
and providing quantitative characterizations (e.g., fluc- 
tuations or full distributions of J's) will be revealing. 
This can also expose the nature of the transition. Does 
it display the same behavior as a typical first-order tran- 
sition? If so, does its hysteresis follow standard proper- 
ties? Does /O c (£ max ,L) scale to /O c (f max /L)? What cor- 
relations, spatial and temporal, can we expect? Struc- 
ture factors, cluster-size distributions, and power spec- 
tra for various quantities should be studied to provide 
us a more complete picture. Moreover, is there a bet- 
ter estimate of the average train length (2H — 1) than 
the rough values from p or (to)? Of course, interesting 
questions also abound for the AC phase (e.g., a system 
with £ max — L). As L — > oo, does Jac attain its max- 
imum at p > 1/2, as our data seem to indicate? If so, 
how can we predict the values of both p and Jac? Sim- 
ilarly, can we understand the significance of the inflec- 
tion point and where is it located? The ultimate goal 
is to find the exact steady state distribution P({rii}), 
likely a highly non-trivial task. Furthermore, we can 
raise all the intricate questions which were directed at 
the ordinary TASEP, e.g., dynamic properties and large 
deviation functionals. Turning to the modeling front, we 
should move beyond the simple AEP and consider more 
complex rules for "pushing" in real applications. Rules 
that readily come to mind include cluster-size dependent 
triggered moves, open AEP, extended particles, inhomo- 
geneous hopping rates, multiple particle species, etc|31]. 
We hope that AEP will open a new chapter for exclusion 
processes as well as a new window into non-equilibrium 
statistical mechanics in general. 
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